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Symmetry and Surface Symmetry Energies in Finite Nuclei 
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A study of properties of the symmetry energy of nuclei is presented based on density functional 
theory. Calculations for finite nuclei are given so that the study includes isospin dependent surface 
symmetry considerations as well as isospin independent surface effects. Calculations are done at 
both zero and non-zero temperature. It is shown that the surface symmetry energy term is the 
most sensitive to the temperature while the bulk energy term is the least sensitive. It is also shown 
that the temperature dependence terms are insensitive to the force used and even more insensitive 
to the existence of neutron skin. Results for a symmetry energy with both volume and surface 
terms are compared with a symmetry energy with only volume terms along the line of /3 stability. 
Differences of several MeV are shown over a good fraction of the total mass range in A. Also given 
are calculations for the bulk, surface and Coulomb terms. 



I. INTRODUCTION 



Moderate to heavy nuclei have neutron excesses due to the growth of the Coulomb energy. The neutron/proton 
ratio N/Z drifts to higher values with increasing A with a value ^ 1.5 near A = 200. The isospin index / = {N — Z)/A 
varies with neutron - proton difference and reaches a value ~ 40/200 — 0.2 near Pb and somewhat higher for the very 
heaviest nuclei. For / 7^ the binding energy of nuclei then acquires an isospin asymmetric part with a dependence 
that is I^A for the volume term. In finite nuclear systems a surface symmetry energy is also present and has a 
^ contribution with dependence I^A"^^^. The energy per particle has a well known expansion in terms of a liquid droplet 
(~| ' model known as the Weizsacker mass formula [l|-|7[. The nuclear symmetry energy along with the Coulomb energy 
determines the nuclear /3 stability line. Nuclei away from the valley of /3 stability will be explored in future rare 
■ isotope or exotic beam accelerator experiments. Systems with much higher neutron excess are encountered only in 



neutron stars which can have values of / approaching unity. In neutron stars only the volume part of the symmetry 
■ energy is important. However, extracting the volume part of the symmetry energy in finite nuclei requires an analysis 
\^ of both the volume symmetry and surface symmetry energy terms. An analysis of the division of the symmetry energy 
00 into volume and surface terms should lead to a better extrapolation of the symmetry energy to the limits involved in 
\l ^ neutron stars. 

i The isospin dependence of binding energies arises from the underlying isospin structure of the nuclear force between 
' nucleons [7| . In particular the force has terms involving the isospin operator f in the isospin conserving form Ti ■ fj 
. which arises from the exchange of isovector mesons. Isospin symmetry is broken basically by the Coulomb force 
B • A small amount of isospin breaking is also present in the nuclear interaction. Nuclear properties explored with 
heavy ion collisions result in systems that are not at temperature T = 0. Non-zero temperatures bring in entropy 
considerations in hadronic systems. Specifically, the Helmholtz free energy F = E — TS becomes an important 
thermodynamic function with natural variables of volume V and temperature T. Nuclear phase transitions are 



' governed by the Helmholtz free energy where energy factors and entropy terms compete to determine the cluster 



distribution seen experimentally |10l . Illj . Isospin properties appear in phase transitions in two component nuclear 
systems. An important example is isospin factionization [10l4l4l | where the dilute gas phase prefers a much larger 
neutron excess than the denser liquid phase. The Coulomb force also plays a role in isospin factionization just as it does 
in determining the nuclear stability line (TB - flTj . Similarly, the phase diagram of nuclei is a function of temperature, 
density and proton fraction which are governed by the interplay of Coulomb and symmetry terms (l5l - [l7| . An extensive 
study of the role of isospin in heavy ion collisions can be found [3 [3 • Discussions of the symmetry and surface 
symmetry energy and how to place limits on the coefficients associated with them can be found in Ref . [20 "22] . The 
role of isospin asymmetry and symmetry energies in nuclear astrophysics are extensively discussed in Ref. [22]. 

This paper is devoted to a study of nuclear energies using density functional theory based on a Skyrme interaction. 
The study involves the following features. A). The nuclear system is finite so that both volume and surface terms 
appear. B). An N ^ Z asymmetry is present so that the more general case of isospin / 7^ is considered with volume 
and surface symmetry energies both present. C). The proton component is charged generating a Coulomb interaction. 
D). The system is allowed to be at temperature T ^ and entropy features are present. The nuclear system is 
therefore non-degenerate. However, the temperature is low enough so that an expansion around the degenerate limit 
can be used. In Sect. 2 the basic relations are summarized. The temperature dependence of nuclear energy of finite 
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nuclei are discussed in Sect. 3 and concluded in Sect. 4. Appendix A gives some semi-analytic expressions for integrals 
used in the density functional theory. 

II. FINITE TEMPERATURE DENSITY FUNCTIONAL THEORY AND THE SYMMETRY AND 

SURFACE SYMMETRY ENERGIES IN NUCLEI. 



A density functional theory will be used to study properties of finite nuclei at non-zero temperature and also 
at temperature T = as a limiting situation. We will limit the temperatures to the low temperature regime so 
that thermodynamic functions can be expanded about the degenerate limit and the T — >■ can be done. The high 
temperature limit is an expansion about an ideal gas and is appropriate for studies of the liquid-gas phase transition 
[isl - flTl . [23I . In this density functional approach a Skyrme Hamiltonian is used which is 
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The H{r) has a bulk part HB{r), a surface part Hs{r) with gradient terms and a Coulomb term Hc{r)- The gradient 
terms are important in finite nuclei and the Coulomb term is important for the charged proton component. The 
^0, til t2, and xq, xi, X2, x^ are parameters. Different choices of these parameters give rise to different Skyrme 
interactions. Here, we consider two Skyrme interactions, SKM(to* = m) and SLy4. These two Skyrme interaction 
have parameter sets given in Table I of Ref. [23| . The m* is the effective mass which is given by 
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At low T or high density, the nearly degenerate proton and neutron Fermi gases have 
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The first term in square bracket is the degenerate limit and the term is the finite temperature correction. The 
Exq is the kinetic energy density, where q = p for protons and q = n for neutrons. 

Since our calculations are done at finite T, entropy S and Helmholtz free energy F become important quantities 
with the connection to the energy E through F = E — TS. The entropy density is 
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for low T or high density to first order in the expansion about the degenerate limit. The density distribution used in 
the evaluations presented in this work is 

Pqc 
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In our evaluation we will first take the proton and neutron radii R to be the same and also the diffuseness parameter a 
to be the same. The central density parameter and Pnc are determined to give the correct number of proton Z and 
neutron N. With this choice of the density we can integrate the energy density of Eq.(IT]) explicitly (see Appendix) 
and thus the energy E(A, Z, T) becomes a function of nuclear size R for a fixed value of diffuseness parameter a. For 
each nucleus with Z protons and N neutrons at a temperature T the nuclear size R can be determined by minimizing 
the Helmholtz free energy. We will also compare this situation with a case where the central densities pqc are the 
same but the proton and neutron radii, Rp and R„, are different. 
The main quantity of interest here will be 

E{A,Z,T) = -B{T)A + Es{T)A^/^ + Sv{T)I^A + Ss{T)I^A^/^ 

where I = {N — Z)/A = {A — 2Z)/A. The E^ij and E^x are the coefficients for the diffuseness correction and 
the exchange correction to the Coulomb energy. For the pairing correction with constant A, c = +1 for odd-odd 
nuclei, for odd-even nuclei, and —1 for even-even nuclei. The above formula at T = is the well known Weizacker 
semiempirical mass formula P, [13, studied extensively by Myers and Swiatecki [26l - l29j . Early studies excluded 
the surface symmetry term Ss and only the surface term Es was included. The values of the coefficients as found in 
textbooks such as Ref.0 are B(0) « 16, Es{0) ~ 17, and Sv{0) « 24 in MeV. The ratio Es/B of surface to bulk 
energy at T = is very close to unity. 

With regard to the volume symmetry term about Y2 the numerical value of the coefficient Sy (T — 0) comes from 
kinetic energy considerations for two degenerate Fermi gases of proton and neutrons. Specifically, the kinetic energy 
contribution is simply related to the Fermi energy Ep as 

SviO)k^n = ^Ef « 12MeV (7) 

The other Y2 of the symmetry energy coefficient arises from interaction terms. For systems with a neutron excess, 
the neutrons and protons experience interaction potentials. In an independent particle model the average potential, 
(called a Lane potential, see Ref.i), felt by a neutron or proton differ. The potential is written in the form 

1 N-Z 

with tz ^ 1/2 for neutrons and ~ ~l/2 for protons. A total potential energy can be obtained by summing the one 
body potential energy and taking Y2 of the result to give 

A depth Vi = 96 MeV gives Sv{0)pot = Vi/8 = 12 MeV, the remaining part of the symmetry energy coefficient 
Sv{0) = 24 MeV. More recently, the importance of a surface symmetry term has been noted. Several calculations 
at r = have been presented regarding this term [20l - [2^ . Considerable variation in the ratio of the T independent 
surface to volume symmetry energy with the difference i?„ — Rp were noted [2^ . Here our focus will be on the 
T-dependent features which arises from the kinetic energy terms and effective mass terms. 

The Coulomb energy Ecoui for a proton distribution of Eq.© (see p. 145 and p. 160 of Ref. 7']) is 



Ec,., = J''''-''c(f, = -e'-[l-{^-') (-) 



5e»(i)"Yv3f^' 



:eV-3— --- — — - (10) 
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with R = roA^/^. For a = 0.53 fm and ro = 1.25 fm, Ec = 0.6912, Edtf = -1.43081, and Eex = -0.5278064 in MeV 
unit. 



III. TEMPERATURE DEPENDENCE OF ENERGY OF FINITE NUCLEI 



To obtain the temperature dependence of the various coefficients, the Helmholtz free energy is minimized for various 
nuclei along (3 stability line with Z protons and N neutrons from ^°Ne to ^°^Pb. The Eq.® is then used to obtain the 
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TABLE I: Energy coefficient minimizing free energy in MeV unit. The rows labeled with "T^ term" are the expansion of the 
coefficient of the explicit T^-dependent term of the kinetic energy coming from the term in Eq.((3l), the kinetic energy term, 
and the rows labeled with "T-indep." are the expansion of the remainders in E{A, Z, T). The rows labeled with "T^{Rn / Rp)" 
are for _R„ 7^ Rp and are explained in text. All other rows are for i?„ = Rp. 





SKM (m* 


= m) 






SLy4 


T (MeV) 





1 


2 


3 





1 


2 


3 


B{T) 


15.310 


15.270 


15.152 


14.955 


15.308 


15.296 


15.264 


15.217 


r-indep. 


15.310 


15.310 


15.310 


15.310 


15.308 


15.308 


15.308 


15.313 


term 




0.039416 


0.039392 


0.039448 




0.011491 


0.011169 


0.010581 


T^iRn 7^ Rp) 




0.039465 


0.039599 


0.039436 




0.011364 


0.011257 


0.010709 


Es{T) 


18.303 


18.804 


20.324 


22.905 


20.008 


20.558 


22.249 


25.205 


T-indep. 


18.303 


18.304 


18.311 


18.340 


20.008 


20.008 


20.024 


20.099 


term 




0.50052 


0.50314 


0.50729 




0.54997 


0.55604 


0.56734 


T^{Rn 7^ Rp) 




0.50037 


0.50240 


0.50738 




0.55043 


0.55572 


0.56694 


Sv{T) 


19.685 


20.108 


21.415 


23.673 


31.113 


31.663 


33.403 


36.656 


T-indep. 


19.685 


19.685 


19.699 


19.745 


31.113 


31.114 


31.152 


31.339 


term 




0.42287 


0.42898 


0.43643 




0.54823 


0.56287 


0.59079 


T\R„ / Rp) 




0.43200 


0.43658 


0.44579 




0.55591 


0.56828 


0.59845 


Ss{T) 


-33.176 


-35.222 


-41.563 


-52.538 


-41.035 


-43.639 


-51.805 


-66.974 


T-indep. 


-33.176 


-33.167 


-33.220 


-33.430 


-41.035 


-41.045 


-41.191 


-41.992 


T^ term 




-2.0551 


-2.0858 


-2.1230 




-2.5948 


-2.6534 


-2.7758 


T^{Rn 7^ Rp) 




-2.0790 


-2.1088 


-2.1460 




-2.5958 


-2.6531 


-2.7905 



expansion coefficients of energy. The results for two different Skyrme interactions, SKM(m* = m) with no effective 
mass and SLy4 with a density dependent effective mass of m* /m = 0.7 at nuclear matter density (23l.[30l|. are presented 
in Tabic U 

The central density pc for SKM(to* = to) has a range 0.14595 - 0.16697 fm'^ at T = 0, 0.14556 ~ 0.16618 fm^^ at 



T= 1 MeV, 0.14446- 



0.16369 fm-3 at T = 2 MeV, 0.14254 

-3 



0.15944 fm-^ at T = 3 MeV, 0.13973 - 



at T = 4 MeV, and 0.13586 - 0.14327 fm"'^ at T = 5 MeV with maximum for f^Ncio 



0.15290 fm-3 

208 T 



The central density pc for SLy4 has a range 0.14286 - 0.14886 fm 
MeV, 0.14017 - 0.14542 fm'^ at T = 2 MeV, 0.13664 - 0.14091 fm 



3 at T = 
3 at T 



and minimum for 8|Pbi26 



0, 0.14223 ~ 0.14799 fm"'^ at T = 
= 3 MeV, 0.13117 - 0.13392 fm^^ 



1 

at 



T = 4 MeV with minimum for ^gf Pbi26 and maximum for ^Zn^^. At T = 5 MeV, it has a range 0.12302 - 0.12519 
with maximum for ^ggSnge and minimum for ^g2Pt)i26 for Z = 30 ~ 82 region and pc — 0.09748 for ipNcio and 
Pc — 0.11736 for |QCa2o- As the temperature T increases, the central density pc which minimizes the free energy F 
decreases and thus the nuclear size R becomes larger. This indicates that no nucleus can be bound with a minimum 
value of F at high temperature due to the entropy which increases as R increases. For ^qNciq with SLy4 parameter, 

there was no pc or R value which minimizes F at T - 5.3 MeV. The density pc which minimizes the free energy F is 
sensitive to the interaction used. 

The results of Table HI from the values of "T-indep." and "T^ term" which are approximately T independent, can 
be summarized approximately at low T as 



in MeV unit for SKM(to* = m) and 



B{T) 

Es{T) 
Sv{T) 
Ss{T) 



B{T) 

Es{T) 
Sv{T) 
Ss{T) 



15.310 -0.0394T^ 



18.303 
19.685^ 



0.50ir^ 
0.423T2 



= -33.176- 2. OSST-' 



(11) 



15.308- 
20.008- 



O.OllST-" 
0.550T2 



= 31.113 -f0.548T^ 
= -41.035 - 2. 595T2 



(12) 



in MeV unit for SLy4. When the Helmholtz free energy F is minimized, the T dependence of energy E(A, Z, T) at 
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low T becomes 

E{A, Z, T) = -(15.31 - QMT'^)A + (18.30 + O.SOT^)^^/^ + (19.69 + 0A2T^)I^A 

-(33.18 + 2.06T2)/2^2/3 ^ Sc-^ + + + cAA-i/2 (13) 

E{A, Z, T) = -(15.31 - 0.nT^)A + (20.01 + 0.55T^)A^/^ + (31.11 + 0.55r2)/2A 

Z^ Z*^'^ 

-(41.04 + 2.59r2)/2A2/3 + i;c^ + + E,,j^ + cA^-i/^ (^4) 

for SKM(m* = to) and SLy4 respectively in MeV unit. 

The temperature dependence of energy Eq.([6]) and Eqs. ljlip - (|14p comes from the T dependent second term in 
Eq.(I31), the kinetic energy. Comparing Eqs.® and @, we can see the entropy TS term is twice that of the T 
dependent part of the kinetic energy with the effective mass E^^^ — ^^^Tq. This means that the T dependence of 

energy originate mainly from the T dependence of the entropy and the T dependence of the free energy F ^ E — TS 
is the same as the T dependence in E with an opposite sign at low T. The entropy can be approximated as 



S = r(o.079A+ 1.00^2/3 + 0.85/2A-4.11/2^2/3\ (^^^^ 



S = Tfo.22yl+ 1.10^2/3 + 1.10/2/1 -5.19/2/12/3) (16) 



for SKM(to* = to) and 



for SLy4 in MeV unit. The volume term of Egs. p^ - (IT51) has the smallest coefficient and the surface symmetry term 
has the largest coefficient among four terms. The minimum free energy F = E — TS of finite nuclei at low T becomes 

F{A, Z, T) = -(15.31 + 0.04T2)A + (18.30 - 0.50T2)A2/3 + (19.69 - 0.42r2)/2A 

Z2 ^2 ^4/3 

-(33.18 - 2.06r2)/2A2/3 + Ec^ + Ed^f^ + E,,j^ + cAA'^/^ (17) 
F{A, Z, T) ^ -(15.31 + 0.11T2)A + (20.01 - 0.55T2)A2/3 + (31.II - 0.55r2)/2^ 

-(41.04 - 2.59r2)/2A2/3 + Ec^ + Ed^f^ + E,,^ + cAA-^'^ (18) 

for SKM(to* = to) and SLy4 respectively in MeV unit. 

The results of Table H] and Eqs. ([TT|) - ([T4l) show the T dependence of the bulk, surface, symmetry and surface 
symmetry energy from T = MeV to T = 3 MeV. At T = the F and E are the same and thus the minimum of 
each of them are also the same. Eq. ^ shows a characteristic T'^ dependence coming from a nearly degenerate Fermi 
gases of protons and neutrons. The bulk and surface terms are similar for the two Skyrme interactions. However, a 
comparison of the symmetry term and surface symmetry term at T = are quite different showing a sensitivity to 



the force used. The large variation of symmetry energy depends on the force used, which is also shown in Ref. 31 1 
at zero temperature. (In Ref. 31] af = Ey/Es is given instead of Es-) The SLy4 interaction has a much larger 
coefficients than the SKM(to* = to) interaction. Most of the difference comes from the T independent part. The 
T2 dependent term of symmetry energy is much less sensitive to the force used. The volume symmetry and surface 
symmetry terms at T = are also sensitive to the difference in neutron - proton radii 122]. However, in Ref.[32j 
an investigation of the neutron - proton radii difference in ^^^Ph using a Skyrme-Hartree-Fock model shows a small 
difference of i?„ — Rp — 0.16 fm compared to nuclear size of about 7 fm. 

The results show that the magnitude of B becomes smaller as T increases while the other coefficients become 
larger. The bulk energy B is insensitive to the temperature having a small coefficient for the T^ term. The surface 
energy Es and the volume symmetry energy Sy have larger temperature dependences which are of the same order. 
The surface symmetry energy Ss is the most sensitive to the temperature having a large coefficient in front of T^. 
The T dependent terms are insensitive to the force used. Most of the dependence on the force used appears in the 
T independent part of energy (see Eq.([T])). This fact shows that the T dependence comes from the kinetic energy 
including an effective mass term, and only the effective mass part has a dependence on force used. The interaction 
part which is sensitive to the force used is insensitive to the temperature. The central density Pc or the nuclear size 
R which minimizes the free energy F is sensitive to the interaction used. From Table HI we can see there is a small 
extra T dependence beside the dominant T^ dependence given by Eqs. dTTI) and ([T^ . Including this small extra T 
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dependence, Ea. pT|) and Ea. p2|) have small corrections that are as follows. By extracting T dependence up to 
directly from B{T), Es{T), Sv{T), and S's(T) in Table U itself, for SLy4, the T dependences of Eq.^ are now 

B{T) = 15.30754 + 0.00122333T- 0.0138600T2 + 0.00114667T3 

Es{T) = 20.00770 + 0.0226783T+0.507110T2 + 0.02095167r3 

SviT) = 31.11260 + 0.0619383T+0.434450T2 + 0.05360167r3 

SsiT) = -41.03514 -0.304690T- 2.059045r2-0.2404950T3 (19) 

The coefficients for SKM(m* — m) have the similar T dependence but less sensitive than for SLy4 case. Specifically, 
Ea.([TT|) is now modified to 

B{T) = 15.30963 - 0.0004933T - 0.0386350T^ - 0.000201667T^ 

Es{T) = 18.30316 + 0.0063483T+0.487620T2 + 0.00720167r^ 

Sv{T) = 19.68525 + 0.0024133T+0.408910T2 + 0.01111667T3 

Ss{T) = -33.17567 -0.011930r- 1.977970T2-0.0564600r3 (20) 

This extra T dependence is due to the different density (central density pc) which minimizes F for different T and 
causes the small variation of the coefficients in Table |T] as T changes and the differences between Egs. dT^ and 
and between Egs. pTj) and (PO)) . The density pc is sensitive to the interaction used as mentioned before. 

Here we used the same size of proton and neutron distribution, Rp = R„. In this case the system has no neutron 
skin. The existence of neutron skin may affects the surface dependence of energy. However since we are interested in 
the temperature dependence of energy we examined the dependent term in kinetic energy (Eq,(|3])) by evaluating 
with different values of Rp and i?„ which are determined by requiring the correct number of Z and N with the same 
central density of ppc — Pnc ~ Pcl"^ where pc is the total central density of the results minimizing F with Rp = Rn- 
The results are shown in Table H] labeled with "T^{Rp ^ i?n)" which are very close to the values of "T^ term". This 
shows that the effect of the different size of neutron and proton distribution is much smaller than the effect of different 
force parameter set. For the case of Rp = i?„ the most asymmetry occurs in the central region while it occurs at the 
surface region for the case of ppc — Pnc- If the overall asymmetry effects of these two extreme cases are the same 
then the total symmetry energy is insensitive to the existence of a neutron skin. The Weizacker expansion of Eq.Q 
cannot distinguish between effects coming from different central densities and different radii. This result is shown 
by comparing the values in the rows of "T^ term" and "T^(i?p 7^ Rn)" for each case in Table |T1 Of course the T 
independent part may have a larger effect for different values of Rp and Rn than the dependent part. A discussion 
of the surface to volume symmetry energy at T = can be found in Ref . ■ 

The energy expansion coefficients for some other Skyrme parameter sets with various effective masses are also shown 
in Table |TT1 In Tables U and |lT] we compared the results for the parameter sets with wide range of effective masses 
from m* /m — 1.0 to 0.577. Eqs.lH)) - ^ show the temperature dependence of energy comes mostly from the effective 
mass. Tables U and [H] show that the above discussions about the qualitative behavior of the temperature dependence 
of energy coefficients are independent to the force parameter sets used. Of course the quantitative values of the T 
dependences depend on the Skyrme parameters used. These tables show that the effect of the effective mass on the 
T dependence of the energy is most visible in the term of the volume energy B{T). For smaller effective mass the 
dependent term of B{T) has the tendency of becoming smaller which can be understood by Eq.([3]). The other 
energy coefficients do not show any specific pattern of the dependence on the effective mass. 

It is of interest to compare the total symmetry energy using the coefficients given in Table U for the two Skyrme 
interactions along the line of /3 stability. In particular we calculate the symmetry energy Egym = Syl^A + SsI^A^^^ 
along a simplified stability line given by 



A 2 1 + 1^2/3- 



(21) 



which is determined by minimizing E{Z, A) = E{A)+asymPA + acZ^/A^^^ with = 0.72 MeV and Usym = 24 MeV 
at T = 0. The = {Np - Zp)/A = {A - 2Z^)/A which is substituted into E,y„^{T) = SviT)PA + Ss{T)PA^/^ 
to give the symmetry energy along the line of /3 stability. The energy is labeled Egymjj — (5y(T)A + SsiT)A'^^'^)I^ 

where Ip = ^^-^ ^"'^ ^ ^ i,^ ^ ^T~a^^) FigUI FigUshows that the T dependence for the symmetry energy is 
more visible for larger A and the symmetry energy is larger for higher T. It shows also that the SKM(m* = m) 
parameter set has a symmetry energy which is too small. 

In Figl2]the difference AEsym,^ between the symmetry energy Esym,i3{T) for SLy4 parameter set and the volume 
symmetry energy with Sy = 24 MeV and ^5 = (dotted curve in FigH]). This figure shows that the symmetry 
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TABLE II: Same as Table |T] but for different Skyrme interaction sets of SkT8 [H, SkM* [s^l, and SkI3 [H with the effective 
mass m*/m of 0.833, 0.79, and 0.577 respectively. 





SkTS 


SkM* 


SkI3 


T (MeV) 


12 


12 


12 


B{T) 
T-indep. 
term 

T^iRr, ^ Rp) 


14.940 14.922 14.869 
14.940 14.940 14.942 
0.01853 0.01810 
0.01851 0.01809 


15.127 15.108 15.051 
15.127 15.127 15.127 
0.01921 0.01911 
0.01918 0.01909 


15.102 15.102 15.104 
15.102 15.102 15.104 
0.00048 -0.00002 
0.00044 -0.00005 


Es{T) 
' T-indep. 
term 

T^{Rn ^ Rp) 


21.684 22.266 24.056 

zi.D04 zi.DoO 2i.<04 

0.5809 0.5881 
0.5810 0.5882 


18.756 19.296 20.947 
18.(56 18. /5/ 18. mO 
0.5388 0.5442 
0.5389 0.5443 


22.386 22.963 24.735 

oo ooc; oo ocy oo AdA 

0.5762 0.5829 
0.5763 0.5831 


SviT) 
1 -maep. 
term 

T^{Rn 7^ Rp) 


24.940 25.506 27.294 

0/1 o/in 0/1 0/10 0/1 OQ1 
z4.y4U .Z4.y4z Z4.yoi 

0.5631 0.5783 
0.5719 0.5871 


29.655 30.210 31.956 

oo fic^p; oo oo (\(^^ 
zy.Doo zy.Do/ zy.Dyi 

0.5527 0.5662 

0.5723 0.5859 


27.382 27.980 29.873 

07 QQO OV QQP; 07 /10V 
Z ( .OoZ Z / .OOO Z / .4z / 

0.5954 0.6115 
0.5982 0.6145 


SsiT) 
T-indep. 
term 

T^{Rn 7^ Rp) 


-20.661 -23.383 -31.920 
-20.661 -20.674 -20.832 
-2.708 -2.772 
-2.717 -2.780 


-43.596 -46.185 -54.325 
-43.596 -43.606 -43.760 
-2.579 -2.639 
-2.627 -2.687 


-27.382 -30.042 -38.438 
-27.382 -27.393 -27.571 
-2.649 -2.717 
-2.643 -2.710 



energy is approximately T independent at around A — 100 with the total symmetry energy of 43 MeV for SLy4 set 
compared to 46 MeV for pure volume symmetry energy only. The symmetry energy for higher T is smaller than the 
symmetry energy for lower T at smaller A than /I = 100 while this result is opposite at larger A. Differences of up 
to 4 MeV were found to be present in the mid mass range of finite nuclei. However for small nuclear masses and for 
a larger nuclear mass {A k, 200 for T = 0) the two choices agree for the total symmetry energy for these coefficients. 
This means SLy4 interaction has the total symmetry energy Egym{T) similar to the value for Esym = 24PA MeV up 
to A^ 200 for T < 3 MeV. 



IV. CONCLUSIONS 

This paper presented an investigation of the energy of finite nuclei using density functional theory. A Skyrme 
approach for the nuclear interaction was used in this study. Bulk, surface, Coulomb and both volume and surface 
symmetry energies are calculated for a nucleus at both zero and non-zero temperatures, ft is shown that the surface 
symmetry energy term is the most sensitive to the temperature while the bulk energy term is the least sensitive. 
Understanding features associated with the volume and surface symmetry energy was the main part of this study. 
Specifically, the volume and surface symmetry energy coefficients Sy and Ss were calculated for two different Skyrme 
interactions, SLy4 and SKM(m* = m), and the results compared. The SLy4 interaction with a density dependent 
effective mass m* is in better agreement with known features of the symmetry energy of finite nuclei. The results 
suggest the importance of a density dependent effective mass. The results also show that the temperature T depen- 
dence of the energy is insensitive to the force used and even more insensitive to the difference of Rp and i?„. These 
results are quite different than the T = results which are sensitive to the force used and the difference of i?„ and 
Rp MM- 

Then a comparison was made between the SLy4 Skyrme interaction results with SviT) and Ss(T) of Table |T] 
{Sv = 31 MeV, Ss = -41 McV at T = 0) and a model without surface effects S's = MeV and with a volume 
coefficient Sy = 24 MeV. The results were compared in Figs[T] and [5] Differences of up to 4 MeV were found to 
be present in the mid mass range of finite nuclei. However for small nuclear masses and for a larger nuclear mass 
{A « 200 for r = 0) the two choices agree for the total symmetry energy for these coefficients. The volume term 
Sv = 31 MeV in SLy4 is much larger than the typical textbook value of Sy ~ 24 MeV of the Weizsacker mass 
formula. The Weizsacker mass formulae with Sy ~ 24 MeV was fit to data without a surface term. Adding a surface 
term that reduces the symmetry energy will result in an increase in the volume coefficient. Various constraints give 
limits [20l-[2^ on Sy with in the range 24MeV < Sy < 33MeV. Understanding the volume coefhcient in finite nuclei 
is necessary in extrapolations of the symmetry energy to neutron star physics. 
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A 

FIG. 1: The total symmetry energy Esym(T) = Sv{T)l'^A + Ss{T)l'^A^'^ in MeV versus mass number A along the line of 
stability. The upper curves arc the interaction SLy4 (the solid curve for T = Q and the dashed curve for T = ,3 McV) and the 
lower curves are SKM(m* = m) (the dash-dotted curve for T = and the dash-dot-dot-dotted curve for T = 3 MeV). The 
dotted curve is a pure volume term symmetry energy with Sv = 24 MeV and Ss = Q MeV. The Eeym,p is in MeV. 



This work supported by Basic Science Research Program through NRF of Korea funded by the Ministry of Educa- 
tion, Sciences and Technology under Grant 2009-0072500 and by Department of Energy under Grant DE-FG02ER- 
40987. 



Appendix A 

For a density p{r) = Pc/{^ + e''""^^/"], using partial integral, 



r'^dr 47r „ r^dr ne^r-R)la 

~ Pc 



[1 + e('-«)/'^]» 2,' " Jq a [1 + e('-«)/'»]"+i 

For the last expression, the lower limit of the integration of y is extended from —R/a to — oo which is a good 
approximation for n > 1 since the symmetric factor e^/(l + e^)^ peaks at y = and becomes zero at ±oo. The 
integral part for various n are 



f 

J — < 





(1 








(1 


+ eyf 







{tt^^R + R^) (A2) 

{-TT^a^'+TT^a^R-SaR^ + R') (A3) 
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A 

FIG. 2: Comparison of a symmetry energy with both volume and surface terms and a symmetry energy with just a volume 
term for SLy4 interaction. The symmetry energy is calculated along the line of /3 stability. The upper and lower figures are 
the difference AEsym,^ and AEgym^p/A in MeV between the two symmetry energies with both volume and surface symmetry 
energy terms for SLy4 interaction and with volume symmetry energy term only (dotted curve in Fi^TJ. The solid curve is at 
T = 0, the dashed curve at T = 1, the dash-dotted curve at T = 2, and the dash-dot-dot-dotted curve at T = 3 MeV. At 
A — 100 and T — 0, the SLy4 interaction gives a total symmetry energy of 43.09 MeV, while the pure volume calculation has 
46.44 MeV. 



dyiay + R?^^^^ = (-fc^^a^ + (6 + n^)a^R - ^ai?^ + r^^ (A5) 

/■°° , , 5ey ( 5(6-|-57r2) . (35 -t- 47r2) , 25 , o\ , 

\^ May + RfjY^ = [- 12 « + 4 « - T"^ + ^ j 

dyiay + Rfj^-^-^^ = (159.37840^ -|- 54.8029a^i? + 7.66445ai?2 + i?3^ (A7) 
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J — oo 



/oo 
dy{ay 
-OO 

/OO 
dy{ay + Rf 
'OO 



(l + e?')5/3 
3j4/3)e^ 

^ {1 + eyy/^ 

(5/3)e^ 

(1 + 62^)8/3 



/■OO 

(8/3)e2' 



/OO 
dyiay + Rf- 
-OO \^ 

/OO 



+ e!')ii/3 



/■°° J / nN3 (13/6)e2' 



(I5.855575a3 + 15.773750^/? + 2.2230575ai?2 + R^) 
(-5.03031250^ + 8.81615625a^jR - 1.335546875ai?2 + R^) 
(-7.80506a^ + 9.10458a^ii - 2.276943aii^ + R^) 
(-11.6424a3 + 10.81947a=^i? - 3.58554oi?2 + R^) 
(-13.26781a^ + 11.836907a^i? - 4.07693333ai?2 + R^) 
(-3.00561a^ + 9.07054a^ii - 0.7352638484ai?^ + R^) 
(-10.7804a3 + 10.331a2i? - 3.30669ai?2 + R^) 



(A8) 
(A9) 
(AlO) 
(All) 
(A12) 
(A13) 
(A14) 



(A15) 



For n = 1, we get the constraint 

f ^3 47r (ay + Rf'cy 4^ . 

J ^ 3 ^V-oo {l + ey f 3 ' 

This relation expresses the central density ppc and pnc in terms of the size R and a given diffuseness parameters a for 
a given number of protons Z and neutrons A'^. Similarly for the gradient dependent term, 

/roo r 1 1 j2 



1 _(_ g{r-R)/a 



dr^ |_l + e(''"-R)/" 



/OO 
dy 
-OC 



(l + e«)2 



(ay + Rf ( e 



1 + ey 

V 



{ay + Rf ( e 



i + ey 



3 2Ea 



1 + 



--2') f- 
3 Mi? 



(A16) 
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